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Abstract 

This paper establishes that generic linear quantum stochastic systems have a pure 
cascade realization of their transfer function, generalizing an earlier result established 
only for the special class of completely passive linear quantum stochastic systems. 
In particular, a cascade realization therefore exists for generic active linear quantum 
stochastic systems that require an external source of quanta to operate. The results 
facilitate a simplified realization of generic linear quantum stochastic systems for 
applications such as coherent feedback control and optical filtering. The key tools 
that are developed are algorithms for symplectic QR and Schur decompositions. It is 
shown that generic real square matrices of even dimension can be transformed into 
a lower 2x2 block triangular form by a symplectic similarity transformation. The 
linear algebraic results herein may be of independent interest for applications beyond 
the problem of transfer function realization for quantum systems. Numerical exam¬ 
ples are included to illustrate the main results. In particular, one example describes 
an equivalent realization of the transfer function of a nondegenerate parametric am¬ 
plifier as the cascade interconnection of two degenerate parametric amplifiers with 
an additional outcoupling mirror. 


1 Introduction 

The class of linear quantum stochastic systems [DEIEIII] represents multiple distinct open 
quantum harmonic oscillators that are coupled linearly to one another and also to exter¬ 
nal Gaussian fields, e.g., coherent laser beams, and whose dynamics can be conveniently 
and completely summarized in the Heisenberg picture of quantum mechanics in terms of a 
quartet of matrices A, B, C, D, analogous to those used in modern control theory for linear 
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systems. As such, they can be viewed as a quantum analogue of classical linear stochastic 
systems and are encountered in practice, for instance, as models for optical parametric 
amplihers [SI Chapters 7 and 10]. However, due to the constraints imposed by quantum 
mechanics, the matrices A, B, C,D in a. linear quantum stochastic system cannot be ar¬ 
bitrary, a restriction not encountered in the classical setting. In fact, as derived in [2] 
for a certain hxed choice of D, it is required that A and B satisfy a certain non-linear 
equality constraint, and B and C satisfy a linear equality constraint. These constraints on 
A, B, C, D are referred to as physical realizability constraints [2]. 

A number of applications of linear quantum stochastic systems have been theoretically 
proposed or experimentally demonstrated in the literature. In particular, they can serve as 
coherent feedback controllers HE], he., feedback controllers that are themselves quantum 
systems. In this context, they have been shown to be theoretically effective for cooling 
of an optomechanical resonator [7], can modify the characteristics of squeezed light pro¬ 
duced experimentally by an optical parametric oscillator (OPO) [H], and, in the setting 
of microwave superconducting circuits, a linear quantum stochastic system in the form of 
a Josephson parametric ampliher (JPA) operated in the linear regime has been experi¬ 
mentally demonstrated to be able to rapidly reshape the dynamics of a superconducting 
electromechanical circuit (EMC) |9]. Linear quantum stochastic systems can also be used 
as optical hlters for various input signals, including non-Gaussian input signals like single 
photon and multi-photon states. As hlters they can be used to modify the wavepacket 
shape of single and multi-photon sources mm- Also, linear quantum stochastic sys¬ 
tems can dissipatively generate Gaussian cluster states ns as an important component of 
continuous-variable one way quantum computers [T3] . 

In certain quantum control problems, such as in coherent feedback H°° [2] and LQG 
[6] control problems, the latter being adapted for addressing an optomechanical cooling 
problem in [7], the important feature of the controller is its transfer function T(s) = 
C{sl — A)~^B + D rather than the system matrices (A, B, C, D). Therefore, an important 
issue in the implementation of coherent feedback controllers is how to realize a controller 
with a certain transfer function from a bin of basic linear quantum (optical) devices. This 
is a special case of the problem of network synthesis of linear quantum stochastic systems 
addressed in [SI El E]. In particular, it was shown in [T3|, generalizing the results of 
DEI E], that the transfer function of all linear quantum stochastic systems which are 
completely passive can be realized by a cascade of one degree of freedom linear quantum 
stochastic systems. Gompletely passive here means that the system can be realized using 
only passive linear optical devices which do not need an external source of quanta for their 
operation. The question of whether cascade realizations exist for general linear quantum 
stochastic systems has remained an open problem. The contribution of this paper is to 
resolve this question by proving that, generically, linear quantum stochastic systems do 
possess a pure cascade realization. This is signihcant from a practical point of view, as it 
allows for a simpler realization of generic linear quantum stochastic systems. 

The remainder of the paper is organized as follows. Section introduces the notation 
and gives an overview of linear quantum stochastic systems and the associated realization 
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theory. Section [^presents a symplectic QR decomposition algorithm. The results of Section 
Inform the basis for a symplectic Schur decomposition algorithm that is presented in Section 
l^and used to show that the transfer function of generic linear quantum stochastic systems 
can be realized by pure cascading. Finally, Section summarizes the contributions of the 
paper. 


2 Preliminaries 

2.1 Notation 

We will use the following notation: i = \/—1, * denotes the adjoint of a linear operator as 
well as the conjugate of a complex number. Vl A= [ajk] then A* = [a*^], and A^ = 
where (■)''^ denotes matrix transposition. 3fJ{R} = [A + A*)/2 and Q'lR} = — A*). 

We denote the identity matrix by I whenever its size can be inferred from context and 
use In to denote an n x n identity matrix. Similarly, O^xn denotes a. m x n matrix with 
zero entries but drop the subscript when its dimension can be determined from context. 
We use diag(Mi, M 2 ,..., M„) to denote a block diagonal matrix with square matrices 
Ml, M 2 ,... ,M„ on its diagonal, and diag„(M) denotes a block diagonal matrix with the 

square matrix M appearing on its diagonal blocks n times. Also, we will let J = ^ ^ 

and Jn = /„ ® J = diag„(J). 


2.2 The class of linear quantum stochastic systems 


Let X = {qi,Pi,q 2 ,P 2 , ■ ■ ■ ,qn,Pn)^ denote a vector of the canonical position and momen¬ 
tum operators of a many degrees of freedom quantum harmonic oscillator satisfying the 
canonical commutation relations (CCR) [qi,Pj] = 2i5ij, [q^qj] = 0) ^ind \pi,Pj] = 0 for 
i,j = l,2,...,n, or, more compactly, xx'^ — {xx'^)^ = 2*J„. A linear quantum stochas¬ 
tic system El Eli G is a quantum system dehned by three parameters: (i) A quadratic 
Hamiltonian H = ^x'^Rx with R = R^ G (ii) a coupling operator L = Kx, where 

A' is an m X 2n complex matrix, and (iii) a unitary m x m scattering matrix S. For 
shorthand, we write G = {S,L,H) or G = {S,Kx, \x^Rx). The time evolution x{t) of x 
in the Heisenberg picture [t > 0) is given by the quantum stochastic differential equation 
(QSDE) (see [HEIE]): 


dx{t) = AQx{t)dt -\- Bq 


;x(0) = X, 


dA{t) 
dA{t)* 

dY{t) = GQx{t)dt-I DQdA{t)^ 

with Ao = 2J„(R + ^{K^K}), Bo = 2zJ„[ -K^S K^S* ], Go = AT, and Dq = S. Here 
Y (t) = (Yi(t),... ,Ym(t))~^ is a vector of continuous-mode bosonic output fields that results 
from the interaction of the quantum harmonic oscillators and the incoming continuous¬ 
mode bosonic quantum helds in the m-dimensional vector A{t). Note that the dynamics 
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of x{t) is linear, and Y(t) depends linearly on x{s), 0 < s < t. We refer to n as the degrees 
of freedom of the system or, more simply, the degree of the system. 

Following |2], it will be convenient to write the dynamics in quadrature form as 

dx{t) = Ax{t)dt + Bdw{t)] a;(0) = x. 

dy{t) = Cx{t)dt + Ddw{t), (1) 

with w{t) = 2(3fJ{A(t)},S{A(t)}, • • • and y{t) = 2(3fJ{Yi(t)}, 

9{Fi(t)},..., ^{Ym(t)}, Here, the real matrices A, B, C, D are in one-to- 

one correspondence with Aq, Bq,Co, Dq. Also, w(t) is taken to be in a vacuum state where 
it satishes the Ito relationship dw{t)dw{t)~^ = (/ + tlrn)dt] see |2]. Note that in this form it 
follows that D is a real unitary symplectic matrix. That is, it is both unitary (i.e., DD^ = 
D = I) and symplectic (a real m x m matrix is symplectic if DSmB>~^ = J^)- However, 
in the most general case, D can be generalized to a symplectic matrix that represents a 
quantum network that includes ideal inhnite bandwidth squeezing devices acting on the 
incoming held w{t) before interacting with the system [U |3]. The matrices A, B, C, D of 
a linear quantum stochastic system cannot be arbitrary and are not independent of one 
another. In fact, for the system to be physically realizable [21EIE], meaning it represents 
a meaningful physical system, they must satisfy the constraints (see HU El El El ID ) 


Mn + + B^rnB'^ = 0, 

(2) 

SnC~^ + = 0, 

(3) 

DJrnD~^ = Jrn- 

(4) 


Note that, more generally, one can consider linear quantum stocastic systems with less 
outputs than inputs by ignoring certain output quadrature pairs in y{t) which are not of 
interest, and a corresponding generalized physical realizability conditions analogous to the 
above can be derived [lU EH]- However, for the purpose of this paper it is sufficient to 
consider systems with the same number of inputs and outputs, as systems with less outputs 
than inputs can then be easily handled [TU] . 

Following [20] 5 we denote a linear quantum stochastic system having an equal number 
of inputs and outputs, and Hamiltonian if, coupling vector L, and scattering matrix S', 
simply as G = (S, L, if) or G = {S,Kx, \x^Rx). We also recall the series produet < for 
open Markov quantum systems [20] dehned by G 2 < Gi — (S 2 Si,L 2 -|- S 2 Ti,ifi -|- if 2 + 
A{L 2 S 2 Ti}), where Gj = {Sj, Lj, Hj) for j = 1,2. Since the series product is associative, 
G„<G„_i<.. .<Gi is unambiguously dehned. The series product corresponds to a cascade 
connection G 2 <Gi where the outputs of Gi are passed as inputs to G 2 ; see [20] for details. 

2.3 Realization theory 

Given an n degree of freedom linear quantum stochastic system with system matrices 
A, B, G, if, how can it be built from a bin of linear quantum components and which com¬ 
ponents are needed? This is the network synthesis question for linear quantum stochastic 
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systems. It was shown in [3] that any linear quantum stochastic system with n degrees 
of freedom system G can be decomposed as the cascade of n one degree of freedom sys¬ 
tem Gi,G 2 ,--- ,Gn together with some bilinear interaction Hamiltonians between them, 
as illustrated in Fig. It was then shown how each one degree of freedom system can be 
realized from a certain bin of linear quantum optical components. 



Figure 1: Cascade realization of Gt, (tt is permutation map of {1,2, ...,?7,} to itself), 
with direct interaction Hamiltonians between sub-systems GttQ) and G' 7 r(fc) for 

j, A; = 1, 2,..., n, following [3]. Illustration is for n > 3. 

In certain control problems, such as H°° and H^/LQG coherent feedback control prob¬ 
lems, it is the transfer function of the systems that is important rather than the sys¬ 
tem matrices G = {A, B,G, D) themselves. The transfer function is dehned as '^g{s) = 
G{sl — A)~^B + D, and rather than realizing a particular quartet {A, B,G, D) one may 
consider realizing Gt = {TAT~^,TB,GT~^,D) for a suitable arbitrary symplectic matrix 
T, since G and Gt have the same transfer function. The transformation T is required to 
be symplectic to ensure that the new internal degrees of freedom z(t) = Tx(t) satishes the 
canonical commutation relations. It was shown in [15] that every completely passive linear 
quantum stochastic system, a system that can be realized using only passive quantum de¬ 
vices, has a pure cascade realization of its transfer function that does require any bilinear 
interaction Hamiltonians between oscillators in the cascade. In the context of Fig. above, 
this means that all bilinear interaction Hamiltonians if 7 r(j) 7 r(j) can be removed. Purely cas¬ 
cade realizations are simpler to implement and are therefore desirable. However, it is not 
known whether the transfer function of general linear quantum stochastic systems outside 
of the completely passive class have such a realization. This is an important open problem 
that is resolved in this paper. 


3 A symplectic QR decomposition algorithm 

The main purpose of this section is to develop a symplectic QR decompositon algorithm 
and derive a necessary and sufficient condition for real square matrices of even dimension 
to possess this decomposition. The symplectic QR decomposition will play an important 
role in proving subsequent results that will be presented in Section]^ We begin by recalling 
some useful dehnitions. 

Let X be an invertible 2n x 2n skew-symmetric matrix and let {o,X«) = o'^X* be a 
skew-symmetric bilinear form on induced by X. A set of linearly independent vectors 
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vi, f 2 ,..., V2n on is said to be a symplectic basis with respect to (o, X*) if (wj, Xvj) = 
— {vj,Xvi) = Xij. The space endowed with (o,X«) forms a symplectic vector space. 
Thus, we shall also refer to (o,X«) as a symplectic form. A matrix T G is said to 

be sympletic with respect to X if XT = X. Therefore, if T is sympletic with respect 
to X, then Tvi,Tv2, ■ ■ ■ ,Tn 2 n is also a symplectic basis for whenever ^ 1 ,^ 2 ,... ,V2n 
is a symplectic basis. In this paper, we will be interested in as a symplectic vector 
space with X = Therefore, unless stated otherwise, it is implicit throughout that 
the symplectic structure on is with respect to the symplectic form (o, J„«). Note the 
standard property that if T is a symplectic matrix then so is and T~^, see, e.g., j21j . 
This property will often be invoked without further comment. 

We also recall the following definition from |15j . 


Definition 1 A square matrix F of even dimension is said to he lower 2x2 block triangular 
if it has a lower block triangular form when partitioned into 2x2 blocks: 


F = 


Fii 

02x2 

02x2 ... 

02x2 


F21 

F22 

02x2 ... 

02x2 

5 

Fnl 

Fn 2 


F 

nn 


ion 

2x2. 

Similarly, 

a matrix F is said to he upper 2x2 


block triangular if F"^ is lower 2x2 block triangular. 


We are now ready to state the main lemma of this section that describes a symplectic 
QR decomposition algorithm with respect to the symplectic form (o, J„«). The lemma is 
based on a symplectic Gram-Schmidt procedure that is different from symplectic Gram- 
Schmidt procedures to construct a canonical symplectic basis in a symplectic vector space, 
e.g., [211 Proposition 40] and its proof. It is in the same class of algorithms as, though not 
identical to, existing symplectic Gram-Schmidt procedures used in numerical analysis with 
respect to the symplectic form {o,]K„«) with ]K„ = J 0 [22]. In fact, our procedure is 

rather analogous to the Gram-Schmidt procedure in spaces with indefinite inner products 
as described in, e.g. [23l Section 3.1]. In the lemma, a symplectic basis is constructed 
sequentially from a given and fixed set of linearly independent initial vectors ui, V 2 , ■ ■ ■, V 2 n, 
which are presented to the procedure sequentially two at a time in that order. As in the 
Gram-Schmidt procedure in indefinite inner product spaces, since the initial vectors are 
given, a certain condition is required for the new procedure proposed below to yield a 
symplectic basis for 


Lemma 2 (Symplectic QR decomposition) Let V be a real invertible 2n x 2n ma¬ 
trix with linearly independent columns vi, V2, ■ ■ ■, V2n from left to right. Let Mi = 
[ V2i-i V2i ] G ]^2nx2 ^ _ 1 ^ 2 ,... , 77 , and Mi = Mi, M2 = [ Ml M2 ], and Mj = 

[ Ml M2 ... Mi] for i = 3, ... , 77 , and assume that Ni = M/JndP* is full rank for 
i = 1,2,... ,77. Then V has a QR decomposition V = SY for some symplectic matrix S 
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and an upper 2 x 2 block triangular matrix Y. Moreover, S can be constructed recursively 
by contructing a seguence of real numbers aj, real 2 (j — 1 ) x 2 matrices Sj, and real in¬ 
vertible 2j X 2j matrices Sj, for j = 2,3, Define pi as the (1,2) element of the 


skew-symmetric matrix Ni, ai = \/\tii\ , 

S\ = 

For j = 2, 3,4,..., n, define Sk recursively as 

/20-1) a 


1 = I 2 , and 

■ 1 0 
0 sgn(/ii) 


S,= [S,., M, ] 


0 




a-i 


1 0 
0 sgn(/ij; 
1 0 
0 sgn(/ij) 


(5) 


T-l 


with aj = \/Jp(j\ and Ej = where fij denotes the (1,2) element of 

{Sj-iEj + Mj)'^^n{Sj-iEj + Mj). Then Sj satisfies SjInSj = Ij for j = 1,2,3,... ,n, 
and S = Sn is symplectic. In particular, the columns of Sk are contained as the first 2k 
columns of Sk+i, thus forming a symplectic basis for k = n. Moreover, defining the 

invertible matrices 




■ 1 

0 

02 x 2 (n-l) 

Xi = 


0 

sgn(pi) _ 


1 

02 (n- 

-l)x 2 

-^ 2 (n-l) 


and 


X,= 


0 

0 


a: 




1 0 
0 sgn(/ij) 
1 0 
0 sgn(/ij) 
0 


0 

0 


-' 2 (n-j) J 

for j = 2,3,...,n, then V = SY for an invertible upper 2x2 block triangular matrix 
Y = X-\ where X = X 1 X 2 ■ ■ ■ X^-iX^. 


Proof. Since is a real 2x2 skew-symmetric matrix and is full rank by hypothesis, 

it is of the form 

= 

with /ii 7 ^ 0. It follows immediately from the given construction of Si in the statement of 
the theorem that S'7Jn>S'i = Ji. Therefore, the columns of Si are mutually skew-orthogonal. 
We proceed further by induction. 

Suppose that the columns of Sj, as constructed according to the theorem, form a 
partial symplectic basis for 1 < j < n, i.e., Sj^inSj = Jj. Consider now the matrix 
Zj^i = SjEjj^i -|- for some real 2j x 2 matrix S 7 + 1 . We will choose Sj+i to satisfy 

Sj 0* This yields the epuation Sj iSnSj Ej+i + SjinMj+i = 0. Since SjinSj = 


0 /ii 

-hi 0 ’ 
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Jj we can solve for Sj+i to obtain Sj+i = Define the real 2x2 skew- 

symmetric matrix Yj^i = We will show that 7 ^ 02 x 2 - We hrst note that 

since Nj+i = [ Mi ... Mj Mj+i Mi ... Mj Mj^i ] is full rank by hypothesis, 

so is [ Sj Mj^i Sj Mj^i ]. This is a consequence of the fact that the columns 

of Sj are, by construction, linearly independent linear combinations of the columns of 
Ml, M 2 ,..., Mj. Moreover, since 

I Si Z,+, 1 = I S, 

while the matrix 

0 I 2 

is evidently invertible, we conclude that 

[ Si ]TJ„| Si 

is full rank skew-symmetric. By the given construction of S'j and S^+i, [ Sj Zj^i Sj Zj^i 

is necessarily of the form 

I s, Z,ii Si Z,+1 ] = 

From the fact that the left hand side of the identity is full rank, it follows immediately 
that Yj_|_i is full rank 2x2 skew-symmetric. Therefore, it is necessarily of the form 

\z — ^ /^i+i 

0 J ’ 

with 7 ^ 0. Dehne cTj+i = ^/\JIj^\ . Consider now the matrix 


Sj 0 

0 K,+1 


0 h ’ 


Zj+I ] 


Zj+i — oijj^iZjj^i 


1 0 
0 sgn(pj+i) 


Some brief calculations shows that, by construction, the matrix Zj+i satishes SjSn Zj+i — 
0 , and i = Ji. Therefore, we have shown for 1 < j < n that if the hypotheses of 

the theorem hold and Sj satishes SjSnSj = Sj then the matrix Sj+i given by (|^ satishes 
5'^iJn‘S'j+i = Sj+i- In particular, S = Sn is a. symplectic matrix. 

Note that by the above construction each Xi as dehned in the lemma is invertible. 
Direct calculations then show that VXi = [ Si V 2 ^n ], VX 1 X 2 = [ S 2 ], • • •, 

VX 1 X 2 ■ ■ ■ Xn = Sn = S, where Vj^n is matrix constructed of columns 2j to 2n of V from 
left to right. Moreover, clearly X = X 1 X 2 ■ ■ ■ Xn is upper 2x2 block triangular since each 
of the Xj in the product has this structure, and therefore so is D = X~^. Hence, V has the 
symplectic QR decomposition V = SX~^ = SY. This concludes the proof of the lemma. 


Let us now look at an example to illustrate Lemma 













Example 3 Consider the matrix 


V = 


-15 

33 

9 

5 


42 

-22 

26 

26 


-12 

7 

-43 

-45 


3 

28 

44 

-37 


It can be verified that V^^ 2 V and the matrix Ni as defined in Lemma\^ are full rank. The 
symplectic matrix produced by executing the symplectic QR decomposition is then 


S = 


-0.4862 

1.0695 

0.2917 

0.1621 


-1.3612 

0.7130 

-0.8427 

-0.8427 


-0.1418 

0.0975 

-0.7193 

-0.7569 


-1.0405 

1.4863 

0.4113 

-1.1093 


while the upper 2x2 block triangular matrix Y such that V = SY is 

■ 30.8545 0 -0.7130 -28.0024 ' 

0 -30.8545 3.2734 -34.7437 

“ 0 0 55.6558 0 

0 0 0 55.6558 


Finally, we give a necessary and sufficient condition for a 2n x 2n matrix V to possess 
a symplectic QR decomposition. 


Theorem 4 Let V and Nj be as defined in Lemma^ Then there exists a symplectic QR 
decomposition V = SY, with S a symplectic matrix and Y an invertible upper 2x2 block 
triangular matrix, if and only if the matrices Ni, N 2 , ■ ■ ■, Nn are full rank. 


Proof. The if part is the content of Lemma For the necessity of the full rankness 
of Ni, N 2 , ..., Nn-i, hrst note that V^SnV = Take any k e 

{1, 2,..., n — 1} and partition F as 


Fll,fc Fi2,A: 

0 Y22,k \ ’ 


with Fii,fc e M2fcx2fc^ Yi2,k e M2fcx2(n-fc)^ Y22,k ^ M2(n-fc)X2(n-fc)^ ^Yieie Fii,fc and Y22,k are 
invertible upper 2x2 block triangular matrices. Since = [ Rk 0 ]F'''JnF[ Rk 0 
from the expression for F above we immediately get that Nj. = F^J^J^Fn ^, which is 
evidently invertible for k = 1,2, ... ,n — 1, while = F'''JI„F is invertible by hypothesis. 


We now provide an example of an instance where the condition of Theorem |^/az/s and 
hence V does not have symplectic QR decompositon. 
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Example 5 Consider the matrix 


V 


0 10 0 

0 0 0 1 
-10 0 0 

0 0-10 


Then V^]inV is full rank. However, it may he verified that the matrix Ni associated with 
V is a zero matrix, hence the condition of Theorem is not satisfied and V does not have 
a symplectic QR decomposition. 


4 Pure cascade realization of the transfer function of 
generic linear quantum stochastic systems 

In this section we employ the results from the preceding section to obtain sufficient con¬ 
ditions under which a (physically realizable) linear quantum stochastic system has a pure 
cascade realization. It is shown that this condition will be met by generic (in a sense 
that will be detailed in this section) linear quantum stochastic systems, thereby extending 
the results obtained in [15] for the special case of completely passive systems to generic 
linear quantum stochastic systems (including a large class of active systems). Let us now 
recall a characterization of linear quantum stochastic systems G that have a pure cascade 
realization, i.e., G can be written as G = Gn < Gn-i < • • • < Gi for some distinct one degree 
of freedom systems Gi, G 2 ,..., G„. 

Theorem 6 115\. Theorem 4] Let R = [Rij]ij=i, 2 ,...,n with Rij G and K = 

[ Ki K 2 ... Kn ] with Ki G C"*. A linear quantum stochastic system G = (S', Kx, Rx) 
with n degrees of freedom is realizable by a pure cascade of n one degree of freedom har¬ 
monic oscillators (without a direct interaction Hamiltonian) if and only if the A matrix is 
similar via a symplectic permutation matrix to a lower 2x2 block triangular matrix. That 
is, there exists a symplectic permutation matrix P such that PAP^ = A, where A is lower 
2x2 block triangular. Let R = PRP^ = [Rij], K = KP^ = [ Ki K 2 • • • Kn ]; with 
Rij G and Kj G C™. If the condition is satisfied then G can be explicitly constructed 
as the cascade connection Gn < Gn-i < ... < Gi with Gi = (S, KiXi, kxjRnXi), and Gk = 

(L, Kj^X}^, 2^^ /br k 2 , . . . , n, where X (^7r(l) ? J^7r(l) ? Q 7 t( 2 ) ^ P-k{2): • • • 5 Qn{n) ? Pn^n )) 7 

and 71 is a permutation of {1,2,... ,n} onto itself such that x = Px. 

Remark 7 The theorem has been stated as a minor and trivial generalization of m 
Theorem 4]- The original did not include the additional freedom of allowing a symplectic 
permutation matrix to transform A into lower 2x2 block triangular form, corresponding to 
a mere permutation of pairs of position and momentum operators in x. For instance, if A 
is in upper 2x2 block triangular form it can be trivially transformed into lower 2x2 block 
triangular form by a suitable symplectic permutation matrix, which by /zs Theorem 4] 
would then be physically realizable by pure cascading. 
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Given an n degree of freedom linear quantum stochastic system G = {A, B, C, D) with 
transfer function Sg(s) = C{sl — A) + D, the problem that will be addressed is how 
to obtain a cascade of n one degree of freedom linear quantum stochastic systems that has 
transfer function Sg(s), if such a cascade exists. Recall that for any symplectic matrix 
T, the system Gt = {TAT~^,TB,GT~^,D) is a physically realizable system that has the 
same transfer function as G. The main strategy is to hnd a symplectic matrix T such that 
Gt is the cascade realization that is sought. Before stating the main result, let us recall 
the real Jordan canonical form of a real matrix; see, e.g., [2H Section 3.4]. Let A be a real 
square matrix then A can always be decomposed as A = VJaV~^, where Ja is a Jordan 
canonical form for A. Of course, although A is real, its eigenvalues and eigenvectors can 
be complex, but they always come in complex conjugate pairs. That is, if A and v are a 
complex eigenvalue and eigenvector of A then so are A* and respectively. Therefore, in 
general, V and Ja may have complex entries. However, when A is real it is also similar to a 
real Jordan canonical form. A real Jordan block for a real Jordan canonical corresponding 
to a real eigenvalue of A is the same as the corresponding block in the Jordan canonical 
form. However, to a pair of conjugate complex eigenvalues A^ = + ibk and X* = ak — ihk 

there will associated with them one or more real Jordan blocks of the upper 2x2 block 
triangular form 

■ Ck 

02x2 
02x2 

. 02x2 

with 


With respect to the real Jordan blocks, A can be written as VJaV~^, where Ja is a real 
Jordan canonical form of A (unique up to permutation of the real Jordan blocks), and V a 
real invertible matrix [241 Section 3.4]. We are now ready to state the main results of this 
section. 


^2 02x2 02x2 • • • 02x2 

Gk I 2 02x2 • • • 02x2 

02x2 Ck 02x2 , 

: ; •. '■ I2 

02x2 02x2 02x2 • • • Gk 


bk 


Lemma 8 (Symplectic Schur decomposition) Let A be a real 2n x 2n matrix. Then 
A has a symplectic Schur decomposition A = S~^US with U lower 2x2 block triangular 
if there exists a real invertible 2n x 2n matrix 

V = [vx V 2 ... V2n-1 V2n ] 


such that 


(i) V brings A into a real Jordan canonical form Ja = V~^AV, with Ja in the upper 2x2 
block triangular form Ja = diag(J^,r, Ja,c), where Jat contains all Jordan blocks cor¬ 
responding to the (possibly repeated) real eigenvalues of A (in upper triangular form), 
and Ja,c contains all real Jordan blocks corresponding to the complex eigenvalues of 

A. 
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The matrices Ni, 

N2, . 

..., Nn-l given by 


Ni = 

[ Vl 

V2 

V2], 


N2 = 

[ 

V2 

V3 Va ]"^Jn[ hi V2 V3 ], 


Nn-l = 

[ Vl 

V2 

••• V 2 n -3 V2n-2V^n[Vl ^2 

■ ■ V 2 n -3 


are all full rank. 

If the above conditions hold, let Si be a 2n x 2n symplectic matrix obtained from 
V by applying the symplectic QR decomposition of Lemma so that V = SiY, 
for an invertible 2n x 2n upper 2x2 block triangular matrix Y as given by the 
lemma, and let P G be a permutation matrix that implements the mapping 

■ ■ ■ ,gn,Pn)"^ (^n, Pn, , Pn- 1 , ■ ■ ■, ?!, Pi )"^ ■ Then A has the sym¬ 

plectic Schur decomposition A = S~^US with S = PSf^ and U lower 2x2 block 
triangular. 

Remark 9 Notice that since A and Ja,c have even dimensions, so does Ja^- One can 
always choose a real Jordan canonical form of A to be of the form Ja = dia.g{JA,ry Ja, c), 
which is upper 2x2 block triangular. 

Proof. Let A = VJaV~^, with Ja and V as given in the lemma. By conditions (i) and 
(ii), using Lemma we can construct a symplectic matrix Si such that V = S'lP, with 
Y real invertible upper 2x2 block triangular as given in the lemma. We can thus write 
A = VJaV~^ = SiYJ aY~^SJ^. Moreover, since Ja = diag(JA,r, </a,c), Y, and Y~^ are all 
upper 2x2 block triangular, the product YJaY~^ is also upper 2x2 block triangular. 
Let P G be the permutation matrix dehned in the lemma. Notice that, by its 

dehnition, the permutation matrix P is symplectic, and that PZP^ is 2 x 2 lower block 
triangular whenever Z is upper 2x2 block triangular. It follows from these observations 
that PYJaY~^P^ is lower 2x2 block triangular. Therefore, we conclude that SAS~^ 
with S = PSj^ is a lower 2x2 block triangular matrix, since SAS~^ = PYJaY~^P^. 
Additionally, notice that S is symplectic since S']”^ (the inverse of a symplectic matrix) and 
P are symplectic. □ ■ 

A direct consequence of Lemma is the existence of the cascade realization of the 
transfer function of a linear quantum stochastic system when the conditions of the lemma 
are satished. 

Theorem 10 Let G = {A, B,C, D) be a physically realizable n degree of freedom linear 
guantum stochastic system. If there exists a matrix V associated to the 2n x 2n matrix A 
satisfying the conditions of Lemma then there exists a symplectic matrix S such that the 
transformed system {SAS~^, SB, CS~^, D) is physically realizable with SAS~^ lower 2x2 
block triangular, i.e., Sg(s) has a pure cascade realization. 
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Proof. Let S be as in Lemma then A = SAS ^ is lower 2x2 block triangular. 
Therefore, 2 ^( 5 ) has a pure cascade realization by Theorem]^ □ ■ 

We emphasize that fulhllment of the full rankness conditions on W, N 2 , ■ ■ ■, Nn-i 
depends on the choice of the matrix V which transforms A into its real Jordan canonical 
form (which is not unique). For some choices of V the full rankness conditions may fail to 
hold and thus a pure cascade realization of the transfer function cannot be obtained. 

Let us call as admissible all real 2n x 2n matrices A satisfying the conditions of Lemma 
and refer to those that do not as non-admissible. The following examples illustrate some 
samples of non-admissible matrices that cannot meet the conditions of Lemma 


Example 11 Consider the matrix 


-1 

0 

0 

1 

_1 

0 

-1 

0 

0 

-1 

0 

-1 

0 

0 

-1 

0 

-1 

T—1 

1 


which has the real Jordan decomposition A = VJaV ^ with (following from Example^ 


- 1 

0 

1 

0 

1 

0 


■ -1 

1 

0 

1 - 

0 

0 

0 

0 

1 


0 

-1 

1 

0 

-1 

0 

0 

0 

j Ja 

0 

0 

1 

1 

0 

0 

-1 

1 

0 


0 

0 

0 

- 1 

T—1 

1 


It can be easily inspected that for any choice of permutation matrix P such that VP satisfies 
the conditions of Theorem one will find that PJaP^ will not be upper 2x2 block 
triangular. This matrix A is therefore non-admissible. 


Example 12 Consider the matrix 

■ -2 0 0 O' 

0-3 0 4 

“ 00 - 10 ’ 
0-40-3 

which has the real Jordan decomposition A = VJaV~^ with 


-1 

0 

1 

0 

1 

0 


1 

1 

0 

0 

1- 

0 

0 

0 

0 

1 


0 

-2 

0 

0 

-1 

0 

0 

0 

j Ja 

0 

0 

-3 

4 

0 

0 

-1 

0 


0 

0 

-4 

-1 

CO 

1 


As with Example 12 it can be verified that fo^r any choice of permutation matrix P such 
that VP satisfies the conditions of Theorem\^ one will find that PJaP^ will not be upper 
2x2 block triangular. Thus A is also non-admissible. 
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We observe the following: 


1. All diagonalizable matrices in 


D 2 nx 2 n 


(including all symmetric matrices) with only 


real eigenvalues are admissible. Example M to be given below involves this type of 
admissible matrix. 


2. Non-admissible matrices in ]R 2 nx 2 n following cases: 


(i) The matrix has a real repeated eigenvalue with geometric multiplicity less than 
its algebraic multiplicity, and there exist two mutually skew-orthogonal real 
basis vectors vi and V 2 for the invariant subspace of associated with that 
eigenvalue. For some matrices of this type it is not possible to permute the 
columns of V and the corresponding rows and columns of Ja to transform them 


into admissible matrices, as illustrated in Example 11 


(ii) The matrix has a pair of conjugate complex eigenvalues A and A* (not necessarily 
repeated) with a corresponding pair of conjugate eigenvectors or generalized 
eigenvectors v and such that the real vectors and —iv+iv"^ are mutually 
skew-orthogonal. Example [T^ illustrates an instance of a non-admissible matrix 
with this property. 


The non-admissibility of a real 2n x 2n matrix entails rather particular properties that 
are unlikely to be possessed by typical matrices. This suggests that admissible real 2n x 2n 
matrices are generic in the set of all real 2n x 2n matrices. Generic is in the sense that 
the set of admissible matrices contains an open and dense subset of 'jpjg jg indeed 

the case and we state it as the following theorem, with the proof being deferred to the 
appendix. 


Theorem 13 The set of admissible 2n x 2n matrices is generic in ]R 2 nx 2 n 


Thus, generic matrices in M 2 nx 2 n pg^yg symplectic Schur decomposition and the trans¬ 
fer function Sg(s) of a generic physically realizable linear quantum stochastic system 
G = {A, B, C, D) has a pure cascade realization that can be explicitly determined us¬ 
ing Lemma and Theorem IT We conclude this section by applying the results obtained 


herein in an example that demonstrates an equivalent realization of the transfer function 
of a nondegenerate optical parametric ampliher (NOPA) by a cascade of two degenerate 
parametric amplihers (DPAs) equipped with an additional transmissive mirror. 


Example 14 Consider a NOPA with two modes Oj = -|- tpj), j = 1,2, satisfying the 

canonical commutation relations [oj, al\ = 5jk and [oj, a^] = 0. The operators describing the 
system is H = f (aja^ — aia 2 ), L = [ y/jai y/^a 2 and S = l 2 - We take 7 = 7.2 x 10’^ 
and e = O .67 = 4.32 x 10'^, values that can be realized in a tabletop optical experiment, see, 
e.g., the experimental work fMSf based on the proposals in The A, B, C, D matrices 
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for the NOPA are: 


-3.6 

0 

2.16 

0 

0 

-3.6 

0 

-2.16 

2.16 

0 

-3.6 

0 

0 

-2.16 

0 

-3.6 


A = 10^ 


B = -8.4853 X 10 ^/ 4 ; C = 8.4853 x 10 ^/ 4 ; D = h 


We can choose the matrix V in Theorem 10 to he 


y = 


0.7071 0 0 0.7071 

0 -0.7071 0.7071 0 

-0.7071 0 0 0.7071 

0 0.7071 0.7071 0 


corresponding to the (real) Jordan canonical form 

= V-^AV = 10^diag(-5.76, -1.44, -5.76, -1.44). 

Using Lemma we compute the symplectic matrix Si and upper 2x2 block triangular 
matrix Y as 


Si = 


0.7071 0 0 -0.7071 

0 0.7071 0.7071 0 

-0.7071 0 0 -0.7071 

0 -0.7071 0.7071 0 

Y = diag(l,-l,l,-l). 


The reguired symplectic transformation matrix from Theorem 10 is S = PSi and 
a cascade realization of the transfer function of the NOPA is Gi = (^i, i?i, Ci, Hi) = 
{SAS-^,SB,CS-^,D) with 


Ai = 10^ 


Bi 

= 10' 

Di 

= h. 


-5.76 0 0 0 

0 -1.44 0 0 

0 0 -5.76 0 

0 0 0 -1.44 

0 -60-6 

6 0 6 0 

-6 0 6 0 

0-606 


; C'l = 10' 


0-6 6 0 
6 0 0 6 

0 - 6-6 0 
6 0 0 -6 


The cascade realization Gi can be decomposed as the cascade Gi = G 12 <^ 11 , with 
Gii = f/2,10' [ ij' “i! 1 [ 1 ,-5.4 X lO^qipi+piqi)] , 
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and 


Gi2 = 



3 3* 

-3 -3z 


^2 

P2 


, -5.4 X 10®(g2P2 +P2g2) 


Each of Gii and G 12 can he realized as a DPA with two transmissive mirrors rather than 
one; see If for details of the realization of Gu and Gi 2 - Note that the pump amplitude for 
each NOPA in the cascade realization is Ax 5.4 x 10® = 2.16 x 10^. Therefore, remarkably, 
the cascade realization G 12 < 6*12 requires less total pump power to realize than the original 
NOPA, i.e., 2 x (2.16 x 10^)^ in the cascade compared to (4.32 x 10'^)^ in the original, 
i.e., half the pump power. So, with the cascade realization one obtains a more power 
efficient realization of the same transfer function which yields the same amount of two¬ 
mode squeezing in the two output beams. 

Finally, note that if V had been chosen differently from the one above, for instance, as 


V 


0.7071 0 

0 0.7071 

-0.7071 0 

0 0.7071 


0.7071 0 

0 -0.7071 

0.7071 0 

0 0.7071 


corresponding to 

= V-^AV = 10^diag(-5.76, -5.76, -1.44, -1.44), 

then it may be readily inspected that the full rankness conditions of Lemma are not 
satisfied, hence this choice ofV cannot lead to a pure cascade realization of the NOPA. 


5 Conclusion 

In this paper we have generalized the ideas and results in mm, that focus on the special 
class of completely passive linear quantum stochastic systems, to show that the transfer 
function of generic linear quantum stochastic systems, which includes a large generic class of 
active systems, can be realized by pure cascading. The proof is constructive as the cascade 
realization, when it exists, can be explicitly computed. This is of practical importance 
as it will allow a simpler realization of a large class of linear quantum stochastic systems 
as, say, coherent feedback controllers or quantum optical filters. Numerical examples have 
been provided to illustrate the results of the paper. In one example, it is shown that the 
transfer function of a nondegenerate optical parametric amplifier has a realization as the 
cascade of two degenerate optical parametric amplifiers having an additional outcoupling 
mirror, which operates for only half of the pump power required by the nondegenerate 
optical parametric ampliher. 

Acknowledgement. Contributions: HN developed the symplectic QR and Schur decom¬ 
position algorithms, associated results and Example SG and IP proved the genericity of 
admissible matrices in discussion with HN. The authors thank the reviewers and Associate 
Editor for their constructive and helpful comments on this paper. 
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Proof of Theorem 13 


Let M 2 n(M) denote the set of 2n x 2n real matrices, M 2 „(M) denote the subset of full rank 
(i.e. invertible) matrices, and M^ 2 n,s(lR) the subset of those matrices in M 2 n{^) that are 
simple (recall that simple matrices are square matrices with simple eigenvalues, i.e., all 
eigenvalues are distinct). M 2 n(M) and M 2 n,s(M) are generic (open and dense) in M 2 n(IR), 
and in fact they are (non-connected) manifolds of dimension 2n x 2n = 4n^. Similarly, let 
A 2 „(M) denote the set of real skew-symmetric 2nx2n matrices, and ^42^(1^) the subset of full 
rank such matrices. A 2 n(l^) is generic (open and dense) in ^42^(1^), and in fact it is a (non- 
connected) manifold of dimension (2n x2n — 2n)/2 = 2n^ — n. For a matrix X in 742n(K), 
we dehne the principal submatrices Xh) as the upper left corner 2i x 2i submatrices of X 
(i.e., the sub-matrices formed by rows and columns 1 to 2i) from up to X^”^ = X. 
Let 742n(K) be the subset of y42n(lR) containing matrices X with the property that all 
Xh),^ = l,... ,77 ,-1 , are full rank. Finally, let M 2 n,s(M) denote the subset of matrices in 
M 2 n,s{^) which are admissible. This means that they have the following properties: (i) 
there is a real invertible 2n x 2n matrix V that puts them in a real Jordan canonical form 
Ja (^ = VJaV~^), which is block-diagonal, with the 1x1 real blocks before the 2x2 
complex blocks (recall that A is simple, so it has no real Jordan blocks of dimension higher 
than two), and (ii) V^]inV G y42n(IR)- The proof uses arguments inspired by the proof of 
genericity of simple matrices in the set of all real square matrices of a given dimension, from 
[28| Section 5.6]. Also, it relies heavily on methods and results from differential topology. 
A standard reference for these methods and results, along with terminology and notation, 
is the book [29]. Finally, a crucial argument uses Theorem 5.16 of |30l Chapter II] and 
its proof. The proof uses Lemma 15 and Proposition and Lemma 16 is needed in the 
proof of Proposition All these results will be proved later on in this appendix. 


Lemma 15 Let T G he a simple matrix with nonzero eigenvalues. There is a neigh¬ 
borhood of T in such that, every matrix in this neighborhood has eigenvectors and 

eigenspaces (the latter represented by projection operators onto the respective eigenspaces) 
which are continuous functions of their entries, and moreover, their eigenvalues are of the 
same type as those ofT. 


Lemma 16 Let F : M 2 n(IR) —>■ A 2 n(lR) be defined by F{V) = V'JnV- Then, F is onto, 
and a submersion (see Section l.f] for terminology). 


Proposition 17 The set ^2n of matrices V G M2n 
open and dense. 


such that V^]inV G A2n 


IS 


We have to show that M 2 n,<j(M) is an open and dense (generic) subset of M 2 n(M). First, 
we show that M2n,s(M) is an open set. Consider a matrix A G M2„,s(]R) C M 2 n,s(M), 
and let A = V.LaV~^. The block-diagonal real Jordan canonical form Ja of A, has the 
1x1 real blocks before the 2x2 complex blocks, and no real Jordan blocks of dimension 


higher than two. Also, V G ^ 2 n- Applying Lemma 15 to A, we conclude that there is a 
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neighborhood N'{A) such that for every A G N'{A), A = VJ^V with V close to V, and 

close to Ja, and with the same block structure. Let E = V — V . Then, V^^inV = {V + 

Eyjn{V + E) = V^JnV + E'^^ry + V^JnE + E'^^nE. Heuce, det = P^KE), a 

multivariate polynomial of degree 4 i in the entries Ej^ oi E = V — V, whose constant term 
is det ■ However, for matrices V G ^2n, det 7^ 0, i = 1,..., n — 1. 

By continuity, there exists £0 > 0 , such that det [V~^SnVY'^'^ 7^ 0 , for any E G M2n(IR) 
with maxjk \Ejk\ < £o- By shrinking the neighborhood of A if necessary, we can satisfy 
maxjfc \Vjk — Vjk\ < So, and hence V G ^ 2 n and A G M 2 „(M). This proves that M 2 n,s(M) is 
an open subset of M 2 n,s(M)- 

To prove that M2n,s(M) is a dense subset of M2n(IR), we must prove that every A G 
M 2 n(M) has a H G M2n,s(M) arbitrarily close to it. Since M2„,s(l^) is a dense subset of 
M 2 n(M), there exists a H G M 2 n,s(IR) arbitrarily close to A. Let A = VJaV~^. Then, Ja 
is a block-diagonal real Jordan canonical form with no Jordan blocks of dimension higher 
than two, and can be structured so that it has the 1x1 real blocks before the 2x2 complex 
blocks. Also, V G M2n(M). If V is not in ^2n, we know from Proposition [Tt] that we can 
hnd a V arbitrarily close to V that is in ^ 2 n- Then, A = VJaV~^ G M 2 n,s(M), and is 
arbitrarily close to A. Hence, M2n,s(M) is a dense subset of M2n,s(l^), and the theorem is 
proven. □ 


Proof of Lemma 15 Theorem 5.16 of 


Chapter H] states that, for a simple ma¬ 
trix in there is a neighborhood of matrices in that contains it, such that the 

eigenvalues of every matrix in this neighborhood are holomorphic functions of the matrix 
entries. Furthermore, in the proof of this theorem, it is shown that the eigenspaces of 
matrices in this neighborhood are also holomorphic functions of their entries. Specializing 
these results to real matrices, we have that for a simple matrix in , there is a neighbor¬ 
hood of matrices in that contains it, such that the eigenvalues and eigenspaces (with 

eigenspaces being represented by projection operators onto the respective eigenspaces) of 
every matrix in this neighborhood are analytic functions of its entries. Let T G 
be a simple matrix with nonzero eigenvalues, and T = U JU~^ a decomposition of it in 
a Jordan form (J is diagonal with distinct entries). Let also T G be a matrix in 

the neighborhood N{T) of T with the aforementioned properties. Taking the entries of 
T arbitrarily close to those of T, the eigenspaces of the two matrices can be made arbi¬ 
trarily close, as well. Hence, we may change the chosen eigenvectors of T (columns of U) 
to form eigenvectors of T in a continuous way. Then, we may write T = UJU~^, where 
U is arbitrarily close to U. Similarly, the diagonal matrix J of eigenvalues of T will be 
arbitrarily close to J. Due to the property of real matrices to have complex eigenvalues in 
conjugate pairs, and the fact that that T has no zero eigenvalues, there is a neighborhood 
N'{T) C N{T) such that every T G N'{T) has eigenvalues not only close, but of the same 
type as T. The reason is that, for a pair of complex conjugate eigenvalues to be created 
(destroyed), two distinct real nonzero eigenvalues must coalesce to a double eigenvalue (be 
produced by the separation of two equal real eigenvalues). This, however, can be prevented 
by shrinking the neighborhood of T as much as necessary. □ 
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Proof of Lemma 16; First, we show that F is properly dehned. Obviously, F{V) is 
antisymmetric. For a P G M 2 „(]R), det P 7 ^ 0, so det F{V) = (det det = (det 7 ^ 
0, and hence FiV) G A 2 n(lR)- Next, we show that F is onto. Let X G y42n(lR)- From 


Subsection 2.5.14], we know that there exists a 2n x 2n orthogonal matrix Q, and a 2n x 2n 
block-diagonal matrix A of the form 


A = diag ^ 


0 Ai 

-Ai 0 


0 An 
-An 0 


with Aj > 0,i = 1,... ,n, such that X = Q"' KQ. Furthermore, the eigenvalues of X are 
±?Ai,..., ±zA„, and hence A* > 0, for a full rank X. It is easy to see that A = D'^^nD, 
with D = diag(-y/Ai, • ■ ■, Then, X = for V = DQ, and V is full 

rank because Q and V are. Hence, F is onto. 

Continuity and differentiability follow from the fact that the entries of F(V) are second 
order multivariate polynomials in the entries of V. Finally, we show that F is a submersion, 
i.e. that its derivative DFy : TvM 2 n(M) Tp(y)^ 2 n(I^) is a surjective linear map from 
the tangent space TvM2n0^) of M 2 n(^) at V, to the tangent space Tp(y)A 2 n(I^) of A 2 n(IR) 
at F(y), see |29l Chapter 1] for terminology and notation. Starting from F(y) = V^]inV 
and “taking differentials” of both sides, we have that dF = {dVY^inV + V^SndV . Hence, 
for a tangent vector v G Ti/M 2 „(M) (infinitesimal variation dV at V), we have DFy{v) = 
v^^inV + V^]inV. The tangent space of M 2 „(]R) at any point V is isomorphic to M 2 „(]R), 
and the tangent space of A 2 n(M) at any point X is isomorphic to A 2 n(lR)- Hence, v G 
M 2 n{^), and DFv{v) G A 2 n(K)- Let tn be a tangent vector in TF(y)A 2 n(IR)- To show 


that DFy : TyM'^ 


2n 


—^ Tp(y\A 


F{V)^2n 


is surjective, we must show that for any such tn, 
there exists at least one v G Ti/M 2 „(M), such that DFy{v) = w. This is equivalent to the 
equation = w having a solution v G M 2 n(M) given any w G A 2 n(K)- Let 

V = in that equation (recall that for any V G M 2 n(M), V is invertible). Then, 

it reduces to —v-\-v'^ = w, where the antisymmetry of Jn, and the identity = —hn were 
used. The general solution of this equation is h = u — ^w, where u is any 2n x 2n symmetric 
matrix. It is to be expected that the solution for v (equivalently for v) is not unique, since 
TvM 2 n{^) is a higher dimensional space from Tp(v')A 2 n(l^)- As a matter of fact, the general 
solution for u, v = (u—is parameterized by a 2 ux 2 n symmetric 

matrix u. The set of such matrices is a linear space of dimension |2n(2n-|-l) = 2n^-|-n, and 
its dimension is exactly the difference of dimensions of TyM 2 nY) (4n^) and TF(y-)A 2 n{ 
(2n^ — n). Hence, we proved that DFy is surjective for every H G M; 

(local) submersion. 


2 n 


I.e. 


F is a 

□ 

Proof of Proposition 17: First, we show that the set of H G M 2 „(M) such that 
F{V) = £ ^ 2 n(I^), is open in M2nY)- Consider such a V. Then, detF(H)(*i 7 ^ 

0,i = l,...,u. Let F G M 2 n(M), and consider detF(I/ + F)^*i. Since F{V + E) = 
{V+Eyin{V+E) = V^InV+E^^nV+V^InE+E^InE, we can see that det FY+EY^ = 
P^YE), a multivariate polynomial of degree Ai in the entries Ejk of F, whose constant term 
is detF(I/)*^®b By continuity, there exists £0 > 0, such that det F(I/ + F)W ^ 0, for any 
F G M 2 n(M) with \Ejk\ < eo,j, fc = 1, 2,..., 2n. This proves that, the set of H G M 2 n{ 
such that F(H) = V^]inV G y 42 n(K), is open in M 2 n 
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Now we shall prove that it is dense as well. It suffices to show that, for V G M 2 n(IR) such 
that F{V) = € (^ 2 n(l^)\^ 2 n(l^)), there exists a G M 2 n(^) arbitrarily close to 

V such that F(V) = V^SnV G ianW- Since X = F{V) = V^SnV G (^„(M)\i 2 n(K)), 
there exists at least one 1 < r < n — 1 , such that detX^'”^ = det F{VY'^'> = 0. Since is 
a 2r X 2r real skew-symmetric matrix, there exists an 2r x 2r orthogonal matrix Q such 
that X^^) = SQ, with 


S = diag 


0 ui 
—ui 0 




5 


and z/j > 0, i = 1,..., r. Then, det = (det QY det S = {ui ... u^Y. Since det = 0, 
this implies that at least one of the z/’s must be equal to zero. Without loss of generality, we 
may assume that the hrst q are equal to zero, 1 < g < r. Let S be given by the expression 
above, where the zero z/’s have been replaced by nonzero Si,... ,eq: 


Let also. 


and 


S = diag 


0 El 
—El 0 

0 Ur 
— Ur 0 



0 

T^q+l 


. -J^q+1 

0 


/ ^ W\^(Q 0 \^( 0 \ 

-W^ U ) \0 h^n-r) ) \ 0 h^n-r) ) ’ 


( 0^ Q \( ~s W\(Q 0 \ 

V 0 hin-r) ) V U )\Q hin-r) ) ' 


It is obvious that X can be arbitrarily close to X, for small enough ei,...,£^, and 
that detX^'”^ = det ^ = (^i ... EqUq+i ... UrY Y 0 - We can also see that detX^ = 
F«(£i,...,eg), for i = 1 ,... ,n, where ... , 5 ^) is a multivariate polynomial of de¬ 

gree at most 2i in the variables ^i,... ,Eq, with constant term equal to detX^). Hence, 
for small enough £ 1 ,... ,Eq, all the determinants detX^) ^ 0 remain so for X. So, by 
slightly changing X to X, we increased the number of principal submatrices X^*^ of full 
rank by (at least) one, compared with those of X. If, for some principal submatrices of 
Xh) (such that detX^*) = 0 ), we still have detX*^*) = 0 , we may apply the same procedure 
sequentially, and end up with a matrix X G y 42 n(K) arbitrarily close to X. From Propo- 
F is globally onto. This guarantees that there exists V G M 2 n(l^) such that 
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sition 

V~^^nV = F(y) = X. Moreover, F is a submersion at V. The Local Submersion Theorem 
m Section 1.4], guarantees that a neighborhood of X = F(y) (in which we may assume 
that X belongs to, because we may construct X to be arbitrarily close to to X) is the image 
under F of a neighborhood of V. Then, there exists V G M 2 n(l^) in said neighborhood of 
V, such that V~^JnV = F(V) = X G y 42 n(IP). Hence, the set of H G M 2 n(IR) such that 
F{V) = G / 42 n(IR), is also dense in M 2 n(M), and the proposition is proven. □ 
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